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Abstract 

Method of invariants is used to obtain effective kp-Hamiltonian with position- 
dependent band parameters and correct boundary conditions for electron and 
hole envelope functions in ^.Bs-lieterostructures with arbitrary interface orienta- 
tion. It is shown that the presence of heterointerface yields additional quadratic- 
in-k and linear-in-k terms in the kinetic energy operator of electron and hole due 
to non-commutativity of position-dependent band parameters and momentum 
operator. In particular, the Tq conduction band Hamiltonian is determined by 
two position-dependent band parameters: a\ = mo/m e and additional parame- 
ter a<i. Similarly for description of the valence band kp-Hamiltonian the three 
Luttinger parameters 71, 72 and 73 are no longer sufficient, and introducing of 
two additional parameters 74 and 75, that determine quadratic-in-k terms, and 
additional parameter fix, determining a linear-in-k term, is necessary. The ad- 
ditional terms in the Tq and T$ kp-Hamiltonians, which are proportional to a<i 
and 74, 75, appear due to interaction at the interface of electron spin and hole 
effective spin with their orbital motion, which is described by envelope functions 
(interface effective spin-orbit interaction). The parameters 02 and 75 arise due 
to relativistic effects and appear in the second order of kp-perturbation theory 
only when spin-orbit splitting of intermediate states is taken into account. The 
parameter 74 is non-relativistic and can be approximately related to the Lut- 
tinger parameters as 74 = (l/3)(— 1 — 71 + 272 + 373). The additional linear-in-k 
term, which is determined by the parameter (3\ , is also non-relativistic and arises 
due to variation of the Bloch functions across the interface. In the framework 
of this approach the r§©r7 valence band kp-Hamiltonian is derived. Presence 
of heterointerface gives rise to short range interface <5-corrections, which are also 
determined by the method of invariants. The new form of the Hamiltonian allows 
us to obtain the correct boundary conditions for envelope functions of electron 
and hole in ^.Bs-heterostructures with arbitrary orientated heterointerface. 

PACS numbers: 71.55.Eq; 73.20.Dx. 

I. INTRODUCTION 
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In spite of the fact that the effective-mass theory for semiconductor heterostruc- 
tures was being widely discussed in the literature [|lj-|§, by now there is no unified 
opinion about the proper form of kp-Hamiltonian and corresponding boundary condi- 
tions for electron and hole envelope functions P|-|[12j|. In the present work it is shown 
that similarly to the case of bulk material the effective mass Hamiltonian in a system 
with heterointerface can be determined in the framework of kp-perturbation theory by 
means of the finite number of band parameters. However, due to absence of a trans- 
lational symmetry these parameters become position-dependent and, generally, their 
number exceeds the number of similar parameters in the bulk material Besides 



that the presence of interface gives rise to additional short range interface ^-corrections, 
which are, in particular, responsible for mixing of light and heavy hole states in the 
T 8 valence band [frjj. The approach used in this paper is close to the one used by 
Karavaev and Tikhodeev || and by Foreman ||, [|TC|j . 

Using the method of invariants we have shown that in contrast to the bulk material 
in a heterointerface system the kinetic energy operator of electron in the T 6 band 
up to quadratic-in-k terms is determined by two position-dependent parameters ati 
and CK2- The first one is the inverse effective mass of electron in the Tq band: ot\ = 
m /m e . The second a 2 is an additional parameter that determines a correction to 
the effective Hamiltonian which is responsible for electron spin-orbit interaction at the 
interface, calculated by means of envelope functions. It arises in the second order of 
kp-perturbation theory when spin-orbit splitting of intermediate states is taken into 
account. Neglecting kp-interaction of the ^(rf) conduction band with all other states, 
except the nearest rg©^^^ valence band, this relativistic parameter is proportional 
to the r^ 5 valence band spin-orbit splitting A so : a>2 = [A so /(3E g + 2A so )](l — m /m e ), 
where E g = - A so = - E^ . 

Similarly it is shown that the hole kinetic energy operator in the T 8 band up to 
the second power of k accuracy (excluding linear relativistic terms) is determined by 
six position-dependent parameters. Three of them 71, 72 and 73 coincide with the 
usual Luttinger parameters for the bulk material, the three additional ones 74, 75 and 
(3\ arise due to the presence of heterointeface. The additional quadratic-in-k terms 
in the T 8 band kp-Hamiltonian, which are proportional to 74 and 75, appear due to 
interaction of hole effective spin (J = 3/2) with it's orbital momentum at the interface, 
calculated by means of smooth envelope functions. The relativistic parameter 75 is 
analogous to the parameter 02 of the Tq band and if we take into account interaction 
of the r 8 (r^ 5 ) valence band only with the nearest r 8 ©r 7 (r^ 5 ) conduction band it is 
proportional to the spin-orbit splitting A' so of the conduction band: 75 = (1/18) (1 + 
71 - 2 72 ) [A'J(E' g + A' so )], where E' g = E^ - A' so = E^ - E^ . The non- 
relativistic parameter 74 can be approximately related to the Luttinger parameters as: 
74 = (1/3) (—1 — 71 + 272 + 373). The parameter (3% determines an additional linear- in-k 
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non-relativistic term and arises due to variation of the Bloch functions across interface. 

The additional quadratic-in-k terms in the T 6 and T 8 kp-Hamiltonians, that are 
proportional to the parameters a<i and 74, 75, which are absent in the bulk material, 
take place only in the non-homogenous systems with disturbed translational symmetry. 
These terms describe interaction of electron spin and hole effective spin with their 
orbital momentum, calculated by means of smooth envelope functions, in the narrow 
interface region. This interface effective spin-orbit interaction takes place in the 
valence band, where the hole should be attributed effective spin J = 1. It explains 
non-relativistic background of the parameter 74 in the T 8 kp-Hamiltonian. 

This approach was used to obtain the hole r 8 ©r 7 (r^ 5 ) two-band kp-Hamiltonian 
in heterosystem. It was shown that if we consider the r 8 ©r 7 valence band states as 
formed exclusively by the r^ 5 state, i.e. neglecting spin-orbit mixing of these states with 
other states, quadratic-in-k terms of the 6x6 hole kp-Hamiltonian are determined by 
five position-dependent parameters 71, 72, 73 and 74. Besides that the T 8 ©r7 valence 
band kp-Hamiltonian includes a non-relativistic linear-in-k term proportional to the 
parameter j3i of the T 8 band. 

Taking into account short range part of the interface potential gives rise to addi- 
tional interface short range ^-corrections in the effective mass Hamiltonian, which also 
include k-dependent terms. These corrections modify boundary conditions for envelope 
functions and their derivatives and lead to additional spin-dependent effects of electron 
and hole interface scattering. Using the method of invariants we can determine the 
form of these corrections for the case of arbitrary orientated hetero interface. 

The effective kp-Hamiltonians obtained in this paper allow us to derive the correct 
boundary conditions for electron and hole envelope functions in /^.Bs-heterostructures 
with arbitrary orientated heterointerface. As an example we obtained boundary con- 
ditions for heterostructure with (001) interface. 

II. EFFECTIVE MASS APPROXIMATION IN HETEROINTERFACE 
SYSTEM. POSITION-DEPENDENT PSEUDO-COORDINATE BASIS 

In contrast to the bulk material, a system with heterointerface does not have a 
translational symmetry and therefore electron wave vector k is not a good quantum 
number any more. The electron and hole scattering processes at the heterointerface, 
changing the wave vector from k to k', lead to the fact that kp-Hamiltonian of the 
heterosystem becomes dependent on two variables k and k'. 

Let's represent the Hamiltonian of a system with heterointerface in the following 
form: 

if = 7r^ + EMx) + /(x).Al/(x). (1) 
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p 2 

1 +C/i(x) 



Here Ui^^x) are the crystal potentials on the left and right sides of the interface 
respectively, which have the same period and symmetry; AC/(x) = £/ 2 (x) — ?7i(x); 
/(x) = /(nx) is a step-like function, which modulates periodic potential in the direc- 
tion perpendicular to the interface; n is a unit vector which is normal to the interface. 
Usually to solve Schrodinger equation with the Hamiltonian (1), the Kohn-Luttinger 
basis \ot,n, k) = e* kx \a,n) is used, where \a,n) is a Bloch state at k = 0, satisfying 
the Schrodinger equation for the bulk crystal with periodic potential t/i(x) 

\a,n) = \a,n) 

2m, 

where a labels bands (including irreducible representations) at the T point of Brillouin 
zone, and n refers to degenerate states (partners). In the Kohn-Luttinger basis the 
Hamiltonian (1) has the following form: 

(a, n, k| H \a', n', V) = (e^ + j£) <W<W<W+ (2) 
/(k - k>)AU™'5 nn , + 4kp^'5 kk , + V^'(k - k'), 1 } 

where /(k — k') = f(k± — k' ± )8k^> is a Fourier transformant of /(x), = nk, 

k|| = k— nk±; p"°,' = (a, n| p \a',n') ; AU aa '5 nn i = AU^(0) are matrix elements which 
determine smooth part of the interface scattering potential and couple states of the 
same symmetry; V^, (k — k') = V^t ^kyk' are interface short range corrections; V^> — 

E f(G ± )AU$(G ± ,G\\ = 0) ; AfO'(G) = / u£(x)Atf(x)u# (x) e ^dx ; G 
Gx^o n 

is the reciprocal lattice vector, G± = nG, Gn = G — nGj_; (1/v / U)m"o( x ) = ( x l Q; ) n ) 
are Bloch wave functions at k = 0, V is the crystal volume, Q is the primitive cell 
volume. 

To apply the effective mass approximation to heterointerface systems it is convenient 
to use pseudo-Wannier basis \a,n,TL m ), which is related to the Kohn-Luttinger basis 
\a, n, k) by a unitary transformation 

\a,n,R m ) = y-]T|a,n,k)e lkRm , (3) 

where R m are discrete Bravais lattice vectors and the summation is performed over 
Brillouin zone. This basis satisfies the conditions of orthonormality and completeness 

(a, n, H m \a , n , R m ') = <5a a '<W<5R m R m ,, 
|a, n, R m ) (a, n, R m \ = 1. 

a,n,R m 
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In the framework of the effective mass approximation it is assumed that the en- 
velope wave functions T/>£(k) = (a, n, k|-?/>) and matrix elements of the Hamiltonian 
(2) (excluding the short-range part) are localized in a small area of k-space inside of 
the Brillouin zone and the summation in (3) can be extended to the whole infinite 
k-space. In this approximation we can consider R m as a continuous variable R and 
the pseudo-Wannier basis transforms into pseudo-coordinate basis 

\a,n, R) = -t=X! I a ' n ' k)e" ikR , (4) 



'V k 

which satisfies the following conditions of orthonormality and completeness: 

(a, n, R|a , n , R') = 5 aa <5 nn i5(R — R'), 

/ dH\a,n,H) (a, n, R| = 1. (5) 

a,n 

The Hamiltonian (1) in this basis according to (2), (4) and (5) has the following form: 



(a, n, R\H\a', n', R') = £ E e ikR (a, n, k\H\a', ri, k') e~ ik R 



£ (a) <W + f(R)AU aa ' 



v k,k' 



6 nn ,6(R - R) - £ [V^5(R - R)] <W<W+ (6) 



4 HV R 5(R - R')] P$ + V^'(R)5(R - R). 



Here /(R) = E /(k)e ikR is a smooth modulating function, V^(R) = EC"'( k )e ,kR 

k k 

are interface short range corrections. In the effective mass approximation we can 
consider /(R) « 0(nR) and ^'(R) « V£?'5(nR), where 6(nR) and 5(nR) are 
one-dimensional step-function and ^-function respectively. 

The short range part of the Hamiltonian, characterized by matrix , is deter- 
mined by the nature of atoms constituting the heterointerface and their disposition. 
For ideal heterointerface the elements of this matrix can be, in principal, calculated by 
microscopic theory |TJ| - |L7|]. In the effective mass approximation symmetry of the 



short range corrections is determined exclusively by the symmetry of the crystal and 
orientation of the interface. Therefore it can be included into the Hamiltonian as an 
additional phenomenological term. Near to the interface plane, which is determined 
by the equation nx = 0, the short range part of the interface potential V s (nx.) can be 
written in the form of the following decomposition: 

V s (nx) = V s '(0) + V s ' (0) Y^niXi + ^ V s " (0) ^ n^jXiXj + ■■■. 

i i,j 
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The position of the interface is determined by the condition V s (0) = 0. Therefore the 
first term in the right hand side of this decomposition can be omitted. The second 
term, which is proportional to V s ' (0), characterizes magnitude of the uniform electric 
field localized in the interface region and is the most significant. The third and the 
following terms characterize non-uniformity of this field. Thus from the phenomeno- 
logical standpoint the matrix of the short range interface correction V can be written 
as a decomposition in the powers of n 

v = Vi n i + ^ E % n i n i h — > 

i i,j 

in which the most significant role plays the term proportional to n. 

The locality of the Hamiltonian (6) allows us, at every point of R-space, to get rid of 
the off-diagonal matrix elements f(R)AU aa ' which mix states of the same symmetry. 
To do that let us introduce a unitary matrix S(R), which at any point R satisfies the 
conditions: S + (R)S(R) = 1 (1 is a unit matrix), and diagonalizes smooth part of the 
potential energy operator of the Hamiltonian (6) 

£ S W '(R) p"") W' + f(R)AU a ' w "] S a " w (R) = E^(R)5 aa , (7) 

a",a'" 

Here E^(R) is position-dependent energy of electron in a band at k = 0, which 
determines the band offset at the interface. Matrix elements S aa (R) are solutions of 
the following system of eigenvalue equations: 

[E (a) 5 aa > + f(R)AU aa '] S a ' a "(R) = E^ a "\R)S aa "(R), 

a' 

where R is a parameter. 

Unitary transformation of the Hamiltonian (6) by the matrix S(R) is equivalent to 
using a new position-dependent pseudo-coordinate basis 

|^,R) = 5»,R> S a ' a (R), 

a' 

which in the effective mass approximation takes account of variation of the Bloch states 
at k = across interface. The Schrodinger equation in this basis can be written as: 

E /rfR / (^,R|^|^,R , )(^',R , |^) = ^(^,R|^), (8) 

a' ,n' 

where the Hamiltonian is 
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S%,R\H\S%,R!)= ]T S +aa {R){a",n,R\H\a'",n',R')S aa {R'). (9) 



Transformation (9) does not break the locality of the Hamiltonian of the heterosystem 
and therefore the equation (8) is a differential equation and can be represented as: 



£ [£^(R)<W<W + I^(k, R)] C'(R) = ^C(R), 



(10) 



where 



(k,R) 



2mo 



m 



+ V c 

1 r n.' 



fRl 



S° 



(R) 



and k = -id/dR, C(R) = (S£,R|V>)- 

Renormalized single-band Hamiltonian H^, (k, R) is obtained by excluding from 
the system of equations (10) all states ^ / a |. Green's function G^'(E,R) = 
1/ E — E^(R) , which is usually used in this procedure, can be approximately sub- 
stituted by G<& [E a {R),R] = 1/ [E a {R) - E^(R)]. Up to the terms of the second 
order this corresponds to the standard perturbation theory. As a result we get 



H^,(k, R) = £W(R)<W + «k, R)+ 



1 rpfid /|* p\ I 

J B( Q )(R)-_B(' 3 )(R)- t mn'l JV ' / 



where omitted right hand side terms correspond to the corrections of the third and 
higher orders of the perturbation theory. 
Commutation relations 



S+(R)k - kS+(R) = i<9S+(R)/<9R, 



S(R)k - kS(R) = i<9S(R)/<9R, 

allow us to transform the single-band kp-Hamiltonian (11) to the canonical form. 
The additional terms, containing <9S + (R)/<9R and <9S(R)/<9R ~ n5(nR), arise due 
to difference of the Bloch wave functions on the both sides of the interface and in 
the framework of the effective mass approximation can be attached to the short range 
(^-corrections. 
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Taking into account all stated above the single-band Hamiltonian (11) can be writ- 
ten as a sum of two terms 

H%(k, R) = H°J«\k, R) + H s n i?\k, R). (12) 

Here (k, R) is a bulk-like part of kp-Hamiltonian which can be written up to the 
terms proportional to kikj in the following form: 



H°^\k,R) = E^(R)6 nn , 



2 m 



h 2 



kp^(R) + p$(R)kj + £r^ M 2 ( R )^> (13) 

J ZTHq ij 



where 



[p!(R) 



Pmn'( R ) 



p^'(R) 



E 



5 +aa "(R)p^'",5 Q '" Q '(R) 



m 0*,, 
m 

is a matrix determining position-dependent band parameters a(R), 



(14) 



(15) 



are position-dependent momentum operator matrix elements, p„°/(R) = p"°/(R). 
In the effective-mass approximation all position-dependent parameters in the bulk- 
like part of the Hamiltonian are proportional to the step-function 0(nR). The sec- 
ond term H^?\k,~R), which vanishes in the bulk material, is a short range part of 
the kp-Hamiltonian. It contains interface short range ^-corrections, including terms 



proportional to rii, riiUj, riikj, riirijki, 



In the effective mass approximation 



<5 2 (nR) « (l/a)5(nR), 5 3 (nR) pa (l/a 2 )5(nR), . . ., where a is a characteristic size 
of the interface region, and all parameters, which determine the short-range part of the 
Hamiltonian, are position-independent. In this approximation 



#2* } (k,R) 



E A^n^nR) + E ^^(nR) + E 



E D^hkjmSinR) + E K^hmSinR)^ + 



+ H.C. 



All these terms can be introduced phenomenologically by the method of invariants, 
based only on the symmetry of the states under consideration (see Section III). 
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To take into account spin-orbit interaction it is necessary to add to the Hamiltonian 
(1) an operator H so , which determines spin-orbit coupling. For heterosystems this 
operator includes three terms which mix states with different spin projections 

H so = 4^ ([VfA(x)p] • <?) + dp/W ([V (AC/(x)) p] ■ a) + 
4^A[/(x) ([V/(x)p] • a) 

The first term in the right hand side of this expression is the usual operator of spin- 
orbit interaction in the bulk crystal with periodic potential t/i(x). The second term 
determines the change of this operator across the interface. The third term is non-zero 
only in a narrow interface region and characterizes spin-dependent interface short range 
corrections. 

Analysis shows that the spin-orbit coupling does not lead to a qualitative change 
of the effective Hamiltonian (12). The dependence E^CR) in (13) still determines the 
interface band offsets, but takes account of the position-dependent spin-orbit splitting. 
The matrix elements of momentum operator p in (13)-(15) must be substituted by 
matrix elements of operator tt, which in the case of a heterointerface system has the 
following form 

All short range corrections, that arise due to the spin-orbit coupling, are relativistic 
and in the framework of the effective mass approximation can be included into the 
short range part of the Hamiltonian (12). 

III. kp-HAMILTONIAN IN HETEROINTERFACE SYSTEM. METHOD 

OF INVARIANTS 

The effective Hamiltonian (12), containing real position-dependent band parameters 
a(R), can be obtained by the method of invariants. To do that it is necessary to perform 
preliminary transformation to k-representation. In this representation the matrix of 
the effective kp-Hamiltonian depends on the three variables k, k' and n, and includes 
Fourier transformants of the band parameters a(k — k') as "constants" of the method 
of invariants. This matrix must satisfy the invariance condition 

E D^Ug^Ug'Xg-^g-^D^ig^ = H^(k,k';n), (16) 

n",n"> 

where D^ a \g) are the T a irreducible representation matrices of the crystal point-group 
F(g e F). The matrix of the Hamiltonian must be also Hermitian 
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^J(k,k»=i/S(k',k;n) (17) 
and besides that it should satisfy the additional conditions of time-inversion invariance 

H%\-k,-V;n) = #2>,k';n). (18) 

Here \a*,n) = T \a,n), where T is the time-inversion operator. 

Since the band parameters a(R) change only in the perpendicular to the interface 
direction, i.e. they are functions only of R± = nR, their Fourier transformants are 

a(k - k') = i / a(R)e- i ( k - k ') R rfR = 
£ /a( J R ± )e-^- fc l)^rf J R ± e- t(k i'- k ii )R i'rfR|| = a(fc ± - tfJV,' 

where R = R — ni?j_. These parameters are invariants under the symmetry transfor- 
mation (16) and also satisfy the additional conditions (17) and (18). 

The method of invariants allows one to obtain the matrix of the effective kp- 
Hamiltonian of the heterointerface system in any order of the perturbation theory 
in k and n with arbitrary "constants" (band parameters) a(k — k'). 

We are mostly interested in A 3 i?5-heterosystems (point-group T d ). Table I. con- 
tains Hermitian combinations of k, k' and n, transforming according to irreducible 
representations of this group, up to the terms proportional kn 2 . 
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TABLE I. Hcrmitian combinations of k, k' and n, transforming according to the 



irreducible representations of the point-group Td- 



Representation 


Even under time inversion" 


Odd under time inversion" 


TV 

1 15 


i(k - k') 


(k + k') 


r, 
1 1 


(k-k') 




Tl2 


< 


2k z k' z k x k' x kyk'y 

{ VS(k X k' x - kyk'y) 




Tl5 




{k y k' z },{k z k' x },{k x k' v \ 




r 25 




i[kxk 


r, , 

1 15 


n 




Ti 


(n • n) = 1 




Tl2 




( 2n 2 z -n 2 -n y 
\ V3{nl-nl) 




Tl5 


^y^z-) ^z^xi *^x'* j y 




Ti 


i{k - k') • n 


(k + k') • n 


Tl2 


< 


I 'I 

{ iVz 


2(k z k' z )n z (k x 
— {ky — k'y)n y 
\_{k x k x )n x {ky 


k' x )n x 
- K)n v ] 


< 


[ 2(fc z 
[ V3[(k 


+ k' z )n z - {k x + k' x )n x 

— {ky + k' y )n y 
x + k' x )n x — {k y + k'y)ny\ 


■p 

J- 15 


1 


' i {{k y - k' y )n z } 
i{{k z - k' z )n x } 
, i{{kx - k' x )n y } 


< 


\ {{ky + k' y )n z } 
{{k z + k' z )n x } 
, {{kx + K)n v } 


r 25 


i [(k - k') x n] 


[(k + k') x n] 


Tl5 




i(k 


- k') • n 2 = i{k - k') 




(k 


+ k') • n 2 = (k + k') 


Tl5 


( i(k x - k' x )(2nl - - ni) 

{ i(K - K)( 2n l ~ n z - n l) 

{ i{k z -k' z ){2nl-nl-nl) 




r {k x 

{ {ky 
{ {kz 


+ k' x ){2nl-n 2 y -ni) 
+ k' y ){2n 2 -nl-nl) 
+ k' z ){2nl-n 2 x -n 2 v ) 


T25 




( i(k x -k' x )(n< -ni) 
\ i(ky - k' y ){n 2 - n 2 x ) 
{ i(kz - K){n 2 x - n 2 y ) 


f {k x + k' x ){n y -n 2 ) 
{ {ky + k' y ){n 2 z - n 2 x ) 
{ {kz + K){nl-nl) 


Ti 


i{k — k') • N 


(k + k') • N 


Tl2 


1 




l{k z -k' z )N z -{k x 

-{ky-k' y )Ny 
\{k X ~ k' x )N X ~ {ky 


- K)N X 

- K)N y ] 


1 


■ 2{k z + k' z )N z - {k x + k' x )N x 

-{ky + k' y )Ny 

^ ^3 [{k X + k' x )N X - {ky + k' y )Ny] 


Tl5 


1 


i{{ky-k' y )N z } 

i{{k z -k' z )N x } 

_ l{{k X -k' x )Ny} 


1 


' {{ky+k' y )N z } 

{{k z + k' z )N x } 

k {{k x + k'jNy} 


r 25 


i [(k - k') x N] 


[(k + k') x N] 



a N x = n y n z , N y = n z n x , N z = n x n y , {AiB^} = \/2{A i B j + AjBi). 



Subsequent transformation to R-representation is performed by the following rules: 
a{k - k')(k ± k') — ► ka(nR) ± a(nR)k, 
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a(k - ^(kik'j) — ► fc i a(nR)jfe j . 
All the terms proportional to rii and n^n, in R-representation contain 5-functions: 

rii — > rij5(nR), 
— >■ n^n^nR), 



M(nR) ± S(nR)ki 



n 



(hi ± k'^rijUi — > kiS(nH) ± 5(nR)fc, 



IV. EFFECTIVE HAMILTONIANS IN A^-HETEROSYSTEMS 



As it was mentioned above a single band Hamiltonian in a heterointerface system 
can be represented as a sum of bulk-like and short-range parts: 

H (a \k, R) = H 0(a \k, R) + # 5(a) (k, R). 

The short range part H s( - a \k, R) contains terms proportional to rij5(nR), 7ij?ij5(nR), 
njfcj<5(nR), rijfcj5(nR)fcj, . . ., which arise in the different orders of perturbation the- 
ory (11). These terms are interface short range corrections which change boundary 
conditions for smooth envelope functions in R-representation. Analysis shows that in 
the short range part of the Hamiltonian it is sufficient to take into account only terms 
proportional to n, n 2 , kn, kn 2 , since these corrections can provide proper form of sub- 
bands dispersion in quantum well, which is dictated by symmetry of the quantum well 
(e.g. point-group D 2 d for a symmetric (001) quantum well). In this approximation the 
short range part of the Hamiltonian can be represented as a sum of four terms 



# s(Q) (k,R) = #5 a) (R) 



H s J a) i 



R) + i/£ a) (k,R) + ^(k,R). 



S(a), 



In this Section the effective kp-Hamiltonians of the r 1; r 6 and T 15 , r 8 bands, and 
also r 8 ©r 7 two-band Hamiltonian are presented. In the bulk-like part H 0( - a \k, R) we 
considered E^(H) = and took into account only terms up to quadratic-in-k. In the 
short range part we included only one term H^ a \H), proportional to the first order 
of n, since it is the strongest short range correction, characterizing a uniform electric 
field localized in the interface region. The rest of the short range corrections are too 
cumbersome and therefore given in Appendix A. 
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It is necessary to note that the bulk-like part of the Hamiltonian if 0(a ) (k, R) in- 
cludes position-dependent parameters, even though their explicit dependence on R is 
omitted. Parameters that describe the short range part of the Hamiltonian H s ^ (k, R) 
are position-independent. 

A. Ti band Hamiltonian 

The I\ conduction band Hamiltonian is determined by one position-dependent 
parameter (inverse effective mass) a = mo/m*: 



re 



2m r 



kak 



E r i-Et, 15 



Here (kak) = k x ak x + k y ak y + k z ak z , a = 1 + ^ 

The short range terms, proportional to the first order of n, in the r x band Hamil- 
tonian do not exist (H^ Tl \R) = 0). 

B. T 6 band Hamiltonian 

The T 6 conduction band Hamiltonian is determined by two position-dependent 
parameters a± and a 2 



re 



2m r 



(ka^k) I + i ( ka 2 k 



kak 



— ki + \aki + 2 — ki + 2dki 



i+l- 



Here J is a unit matrix 2x2,^ are Pauli matrices, 

The parameter ol\ = m /m* determines the electron inverse effective mass in the T e 
band and is equal to the parameter a of the Ti band, when we neglect spin-orbit cou- 
pling. The parameter a 2 determines an additional relativistic term, which arises when 
the presence of hetero interface is taken into account. It vanishes in the bulk material 
due to commutativity of a 2 and operator k. This term describes electron effective 
spin-orbit coupling at the heterointerface, which is calculated by means of envelope 
functions (interface effective spin-orbit interaction). Indeed, using the commutation 
relation a 2 (R)k — ka 2 (R) = 2<9a 2 (R)/<9R and taking into account the fact that in the 
effective mass approximation <9a 2 (R)/<9R A« 2 n5(nR), this term can be represented 



in the form A« 9 



n x k 



er) 5(nR). It is easy to see that it is analogous to the op- 
erator of spin-orbit interaction, in which the factor proportional to Aa 2 no~(nR) acts 
as an effective electric field localized in the interface region and k is the momentum 
operator in the effective mass approximation. 

The parameter a 2 can be approximately determined if we take into account kp- 
interaction of the r 6 (I^) conduction band only with the r 8 ©r 7 (r^ 5 ) valence band. In 
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this approximation it is proportional to the spin-orbit splitting of the T\ 5 band and to 
the square of the matrix element (r^ 5 ||p|| 1^), which determines the electron effective 
mass in the r 6 (r^) band. Thus it can be approximately represented as 

a 2 = [A so /(3E g + 2A so )] (1 - m /m e ) , 

where E g = E^ 6 — E^ 8 , A so = E^ 8 — E^ 7 is the spin-orbit splitting of the r^ 5 band. 

The short range part, proportional to the first order of n, in the band Hamilto- 
nian does not exist (iff (re) (R) = 0). 

C. Ti5 band Hamiltonian 

The bulk-like part of the r 15 band hole Hamiltonian, containing four quadratic-in-k 
terms and one linear-in-k term, can be written in the form 



(k&xk) / - 6 E [khk ~ | (k6 2 k) 



12E{fci6 3 fci+i}{^+i} + *3( 



k6 4 k 



•J) + -^V3E 



tt(> 



i 2 - 

{Ji+lJi+2} 



where 7 is a unit matrix 3x3, Jj are matrices of the angular momentum J = 1, 



= kid — dki, ao is the 



{JiJj} — 2 ( t ^i t ^7 JjJi)i ^k%oikj\ — 2 (j^iClkj -\- kjQ,kj^j , k{d 

lattice constant, which was introduced for convenience. 

Band parameters bi can be expressed by position-dependent parameters p/, which 
describe kp-interaction of the Ti 5 valence band with the Ti, Ti 2 , and T 2 5 bands: 



61 = -1 + 2pi + 4p 2 + 4p 3 + 4p 4 , 



&2 = Pi - 92 - P3 + 2p 4 , 



h = Pi - 92 + P3 - P4, 



b 4 = pi + p 2 - p 3 - Pa- 



Here 



1 >P K£i5 ||p|| 
Pl 3m V E^-E^s ' 
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1 K r i5 IIpII r 25;g)| 2 



6m Y E g 25 - 



i K r i5 IIpII r i5;g)| 2 



6m Y E q ls ~ EVk 



= 7— S 



I(ri5 IIpII r 12 ; g>| 2 



6m Y E [ 12 - EFl5 

The position-dependent parameters &i, hi and 63 are analogous to the Ti 5 band 
parameters of the bulk material. The parameter 6 4 arises due to the presence of the 
heterointerface and similarly to the parameter a<i of the T 6 band describes interface 
effective spin-orbit interaction. In this case hole in the T 15 band should be attributed 
an effective spin J = 1. If kp-interaction of the valence band with all the T25 bands 
is neglected (p2 = 0) the parameter 64 can be expressed by all other parameters b\, b 2 
and 63. In this approximation 64 = (1/3) ( — 1 — b\ + 26 2 + 363). 

The additional linear-in-k term, which is proportional to a position-dependent pa- 
rameter d, appears in the T 15 band kp-Hamiltonian only when the presence of interface 
is taken into account and arises due to variation of the Bloch functions across the inter- 
face. Indeed, the diagonal matrix elements of the momentum operator p^ 1 " 15 = pi^ 1 / 5 " 1 
are equal to zero since the operator p is odd under time inversion. However, this does 
not take place when the position-dependent pseudo-coordinate basis is used, since in 
this case, according to (15), the diagonal position-dependent matrix elements of the 
momentum operator also include off-diagonal matrix elements of the operator p, which 
are non-zero in general 

p!S 5) (R)= E 5+ r - r '-(R) P S r -^'^-(R). 

TV TV/ 

1 15 1 15 

The short range part of the Tis band Hamiltonian, which is proportional to the first 
order of n, is determined by one position-independent parameter r and has the form 

h 2 

H^\R) = 2v^r5>{<Wi +2 }5(nR). 

zm a ■ 



This short range correction coincides with the correction, introduced in Ref. flljj. for 
the case of (001) interface, and generalizes it for arbitrary interface orientation. 

D. T 8 band Hamiltonian 



15 



The bulk-like part of the r§ band hole Hamiltonian along with the usual quadratic- 
in-k terms, determined by position-dependent Luttinger parameters 71, 72 and 73, is 
also characterized by two additional quadratic-in-k terms, determined by parameters 
74 and 75, and two linear-in-k terms, described by parameters fa and fa 



^)(k,R) = ^ 



k 7l k) / - 2 E \kl2k - \ (k 72 k)l Jf- 



4 E {klzk+x} Wi+i} + «2 ([k 74 kl • J) + 28 E [k 75 kl Jf+ 

i^E [hfa] {Vwl + i^{kfa} {J? +1 - Jf +2 )} 



x 



4, Jj are matrix of the angular momentum J = 3/2, 



the Tg band parameters can be determined by 
7i = h, 72 = b 2 , 73 = h, 74 = &4- The addi- 



where J is a unit matrix 4 
jfcjaj = I (/^a + ak^j. 

Neglecting spin-orbit interaction 
the T 15 valence band parameters bf. 
tional parameters 74 and 75, which describe quadratic-in-k terms, arise due to interface 
effective spin-orbit interaction, calculated by means of envelope functions (the hole ef- 
fective spin J = 3/2). The parameter 74 is non-relativistic and can be determined 
approximately by the Luttinger parameters if kp-interaction of the r^ 5 band with all 
the T 2 5 bands is neglected (i.e. p 2 = 0). In this case 74 = (1/3) ( — 1 — 71 + 272 + 373). 
The first four terms of the T 8 valence band kp-Hamiltonian, which are proportional 
to the parameters 71, 72, 73 and 74, coincide with kp-Hamiltonian obtained by Fore- 
man in Ref. 
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The parameter 75 is non-zero only if spin-orbit splitting of the 
intermediate states in the kp-perturbation theory is taken into account. This param- 
eter can be approximately determined, by taking into account kp-interaction of the 
r 8 (r^ 5 ) valence band only with the nearest r 7 ©r 8 (r < j; 5 ) conduction band. Then it is 



7s 



;i/18)(l + 7 i- 2 72 ) A'J(E'+A' so 



where E' g 



E r \ A' 



is the spin-orbit splitting of the r^ 5 conduction band. 

The hole Tg band Hamiltonian contains two linear-in-k terms, determined by the 
parameters fa and fa. The parameter fa, which arises due to variation of the Bloch 
functions across the interface, is non-relativistic and is originated from the correspond- 
ing parameter d of the T 15 band (fa = d if spin-orbit coupling is neglected). The 
parameter fa is relativistic and analogous to the corresponding parameter of the r§ 
band in the bulk material, that arises due to the absence of the inversion symmetry in 
y4 3 5 5 -semiconductors. 

The short range part of the T 8 band Hamiltonian, which is proportional to the first 
order of n, is determined by one parameter u, which is originated from the correspond- 
ing parameter r of the Tis band (when the spin-orbit coupling is neglected u = r) and 
has the form 
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E. Two-band r 8 ©r 7 Hamiltonian 

The r 8 ©r7 two-band hole Hamiltonian can be represented in the following matrix 
form: 



iJ( r ^) + (k,R) H^(k,R) 

Parameters that determine this Hamiltonian can be approximately related to corre- 
sponding parameters of the r 8 (T^ 5 ) band Hamiltonian if we neglect spin-orbit mixing. 
In this case H^ Ts \k, R) coincides with the T 8 band Hamiltonian (see Section IV. D), in 
which 75 = f3 2 = 0. Note that the T 7 band Hamiltonian also includes a term which de- 
scribe interface effective spin-orbit interaction. However in contrast to the r 6 (r^) band 
Hamiltonian this term is non-relativistic and proportional to the parameter 74. Also in 
this approximation the H 0( - rsr? \k, R) block contains one linear-in-k term proportional 
to the parameter j3\ 



H°^\k,R) = A so I + 



-.2 r 



2m 



(k 7l k)/ + i2([k7 4 k 



a 



VS - k yl2 k y ) J 2 ri2 ] {ki+ilsh+2 } If 15 + 

% E [^ + i74^ +2 ] I? 5 + ^75 E [hPi] I? 

where / is a unit matrix 2x2, A so = E^ 8 — E^ 7 , matrices if" are given in Appendix 
B. 

The short range part of this Hamiltonian, which is proportional to the first order 
of n, can be described by the corresponding parameters of the r 8 band 

ff„ s(r '>(R)=0, 
V. BOUNDARY CONDITIONS 
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The boundary conditions for envelope functions can be obtained by integrating the 
effective-mass equation 

E^2(k,R)#(R)=^Vi a) (R), 

n' 

along the normal n across the interface. Taking into account in the effective-mass 
Hamiltonian H^(k,R) only bulk-like part H 0( - a \k,R) and main short range term 
H^(\i, R), proportional to the first order of n, we obtain the following boundary 
conditions, that depend on the interface orientation 

V'l nR=+0 = ^1 nR=-0, 
(AiR^O |nR=+0 - (A n Rll>) |„R=-0 = ^nR^ | nR=0 , 

here matrix operator A nR is determined by the bulk-like part H 0( - a ^ (k, R) of the Hamil- 
tonian, and matrix B nR by its short range part i?^ a ^(k, R). 

The matrices A nR and B nR , that determine boundary conditions for all the bands 
discussed in this paper, are given below for the case of (001) interface, when n x = n y = 
0, n z — 1 (nR = z). 

A. Yi band boundary conditions 



B< r >> = 0. 
B. T 6 band boundary conditions 



M 6) = 



h 2 



2m 



(a-Jz^ I + ia 2 (-(J x k y + a y k x ^ 
B { ? 6) = 0. 



C. Ti5 band boundary conditions 



18 



4(ri B ) = 

z 2m 



(b 1 k z )l~2b 2 (2JI-JI-JI) k- 
Qb 3 {{JyJz} ky + {J Z J X } k^j + i3b 4 [~J x k v + J y k x ^ + ^V^d {J x J y } 



D. T 8 band boundary conditions 



i(r 8 ) = J±_ 

z 2m 



(71**) ' " §72 (2J 2 2 - J 2 - JJ) fc,- 



273 



j__2 



E. r 8 ©r 7 band boundary conditions 



^(r 8 er 7 ) 



_g(r 8 er 7 ) 



i(r 8 ) i(r 8 r 7) 
i(r 8 r 7) + ^(r 7 ) 

5(r 8 ) #(r 8 r 7 ) 
^(r 8 r 7 )+ ^(r 7 ) 



where A^ 8 ) and coincide with corresponding matrices of the T 8 band (see Section 
V. D.), in which 75 = and /3 2 = 0. 



h 2 



2m n 



^(r 8 r 7 ) = ?£_ 

z 2m 



(71^) I + 2274 (-(T x k y + cr y k x ^ 

-2V2 l2 I^k z - VQls (i^'ky + I^'h) + 
^{-^h + iyk x ) + ij= 2 PiI* 15 
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£f 7) = 0, 



2m a y/2 z 

The k-dependent short range terms, which were not included into the boundary condi- 
tions, along with ^-functions contain also their derivatives, which modify the boundary 
conditions. However, influence of the omitted corrections is small and therefore when 
considering systems with heterointerface (e.g. quantum well) it is more convenient to 
take these corrections into account by the perturbation theory in the framework of the 
envelope function approximation. 



V. SUMMARY 



In this paper using a position-dependent pseudo-coordinate basis the effective mass 
equation for electrons and holes in heterosystems for arbitrary interface orientation 
is obtained. Using this basis one can exclude off-diagonal matrix elements of the 
interface scattering potential and in the explicit form introduce into the theory position- 
dependent parameters, that describe bulk material (electron energy at k = and band 
parameters) and can be taken from experiment. This approach is close to the theories 
developed by Karavaev and Tikhodeev B and Foreman M, WU . 



To obtain the effective kp-Hamiltonians with position-dependent parameters in 
74355-heterointerface systems a modified method of invariants is suggested. Since the 
presence of interface disturbs translational symmetry of the crystal and leads to scatter- 
ing of the particles, changing their wave vector from k to k', the effective Hamiltonian 
of the heterosystem, in contrast to the bulk material, should be considered dependent 
on these two variables. Besides that the effective kp-Hamiltonian has to depend on 
the interface orientation, and therefore it should be described by one more indepen- 
dent variable n (the normal to the interface). Using this approach we obtained for 
arbitrary interface orientation the effective kp-Hamiltonians with position-dependent 
band parameters for the Ti and Tq conduction bands, and T$ valence bands, and 
also r 8 ©r 7 two-band effective kp-Hamiltonian. 

In the framework of this approach the effective kp-Hamiltonian with position- 
dependent parameters can be represented as a sum of two parts: bulk-like part if 0(a ) 
and short range part H s ( a \ The bulk-like part, along with the usual terms, that 
are present in the Hamiltonian of a homogenous material, also includes additional 
quadratic-in-k terms. They arise due to non-commutativity of the momentum opera- 
tor and band parameters and vanish in the bulk material. Analysis shows that these 
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terms describe effective spin-orbit coupling at the interface and lead to spin-dependent 
effects of electron and hole interface scattering. 

Besides that the bulk-like part of the kp-Hamiltonian iJ 0(Q ) can include non- 
relativistic linear-in-k term, that appears due to variation of the Bloch functions across 
the interface. The short-range part H s ^ includes terms proportional to n, n 2 , nk, nk 2 , 
. . .. The terms proportional to the first order of n, which characterize the uniform elec- 
tric field localized in the interface region, are the most important. In particular, in the 
Tg band these corrections lead to a well known effect of heavy-light holes mixing at the 
interface. Analysis shows that kn and kn 2 type terms are responsible for linear-in-k 
Rashba-like terms in the electron and hole subbands dispersion in quantum wells. In- 
fluence of these terms on the effect can exceed the contribution of the cubic-in-k terms, 
which exist in the bulk Hamiltonian of the A 3 i?5-semiconductors due to absence of the 
inversion symmetry. 

The effective Hamiltonians derived in this paper allow one to obtain the correct 
boundary conditions for electron and hole envelope functions for arbitrary interface 
orientation. The additional terms, which were introduced into the Hamiltonians in 
this paper, modify usually used boundary conditions. However, some short range 
corrections (for example n 2 , nk, nk 2 type terms) are small and can be easily taken 
into account by the perturbation theory in the framework of the envelope function 
approximation. 



APPENDIX A: HIGH ORDERS SHORT RANGE PARTS OF THE 
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HAMILTONIANS 



Al. Ti band Hamiltonian 




2m ao 



#(n-n)5(nR), 





Here IkS] = k5(nR) - 5(nR)k. 
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A2. T 6 band Hamiltonian 



H S k ^\k,R) = 



H S f*\R) = 
h 2 



-f|(n-n)/5(nR), 



2m n 



2m Q a 

iAi([k<j] ■n)/ + A 2 ([{k5} xn 



a 



2mo 



([ktf] • N) I + n 2 E {US} (n 2 +1 - n 2 i+2 )a t + 
^([{wJxNj.ff) 

Here {k5} = ± [k5(nR) + 5(nR)k], / is a unit matrix 2 x 2, a% are Pauli matrices. If 
we neglect spin-orbit mixing rj — g, X± — I, /ii — q and A2 = /X2 = ^3 = 0. 

A3. Ti5 band Hamiltonian 



si(n • n)/5(nR) + 6s 2 EK? - \)J 2 5(nR)+ 



TT°( L 15j/D\ 

12s 3 EiVa<W* + 2}5(nR) 



^i ri5) (k,R) = £- 



2mo 



it! ([ktf] • n) / + i6t 2 E ([M] n< - § ( 
2l2t 3 E { [M] n m } {J^+i} + 3t 4 ([{k5} x n] • j) 



n U?+ 



^ 15) (k,R) 



2mo 



i2A/3wi E M (n • n) { J i+1 J i+2 } + 



i2VSw 2 E [M] (2n? - n 2 i+1 - n 2 +2 ) {J l+ iJ t+2 } + 
3u> 3 E {M} - n 2 +2 ) Ji + iw 4 ([ktf] • N) J+ 
^6w 5 E ([kiS] N t - § ([k(j] • N)) Jf+ 

H2w 6 E{[ks] N i+1 } {JiJ i+1 } + 3w 7 ([{k5} x N] • J) 

Here / is a unit matrix 3 x 3, Jj are matrices of the angular momentum J = 1. 
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A4. T 8 band Hamiltonian 



R) = 



2moao 



6(n ■ n)/5(nR) + 2£ 2 £(ra? - |) J?5(nR) + 



46EiVa^+i^ + 2}5(nR) 



^T 8) fk,R) 



a 2 

2mo 



[ktf] •n)/ + i2r 2 E([M 



i4r 3 E {[M] ™m} {-Vi+i} + 2r 4 ([{ktf} x n 



k5 

■jH 



n U 



8r 5 E [{&} x n]jf + ^r 6 E {{ W} n i+2 ) {j t - J? +2 ) } 



#££ 8) (k, R) = ^ i^xi E M ( n ■ n ) + 

275X2 E [4^] (2n? - n? +1 - n? +2 ) { J m J i+2 } + 
2X3 E {ks} (n? +1 - n? +2 ) Jj + % ([k5] • n) /+ 
i2* 5 E ([M] k - \ {[kS] ■ N)) J 2 + iA X6 E { [ks] N i+1 } {J t J t+1 } + 
2 X 7 ([{k5} x N] • J) + -^ X s E (n ■ n) { J, (j 2 ^ - J^ 2 )} + 

^X9 E {45} (2nf - n 2 +1 - n? +2 ) { J, ( Jf +1 - Jj +2 ) } + 8x10 E {U} (<i - < 2 ) J?+ 

^Xn E { {fc+i*} (^1 - JI2) } + 8X12 E [{k<f } x n] . Jf . 

Here / is a unit matrix 4 x 4, Jj are matrices of the angular momentum J = 3/2. 
If we neglect spin-orbit mixing £j = s$ (i = 1,2,3), Tj = U {i = 1,2,3,4), — w % 
(i = 1, . . . , 7), r 5 = r 6 = xs = • • • = X12 = 0. 

A5. r 8 ©r 7 band Hamiltonian 



H, 



5(r 8 er 7 ) 



(k,R) 



#f (rs) (k,R) 



H, 



5(r 8 r 7 )+ 



i/f rsr7) (k,R) 

s<r 7 ), 



(k,R) tf|^(k,R) 



Here £ = n 2 , /cn, /cn 2 ; iff^(k, R) coincide with the corresponding matrices of the T 8 
band (see Appendix A4), in which r 5 = r 6 = xs = • • • = X12 = 0. 
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^T 7) (k,R 



#2 r7) (R) 
h 2 



-ei(n-n)/5(nR), 



2m r 



^ 7) (k,R) 



2m a 

IT! ( k5 ■ n) I + 2r 4 ( {k#} x n • a) 
2X3 E{M} (^ +1 -nf +2 )a,+ 



2mo 



f7-5(r 8 r 7 ) 



R) 



%([y].N)/ + 2 XT ([{k(S}xN]. ( 7) 



2 mo ao 



v^EA^^nR) 



if; 



s(r s r 7 ) 

kn 



(k,R) 



n 2 

2mo 



iV2r 2 [(2 



n 2 - 



>, x - [M] n J/ )/ 1 ri2 + v / 3([U" 
n, +2 } if 15 + ^r 4 E [{k*} x n^/f 25 



^5 



n,l/ 2 ri2 



+ 



H 



s(r 8 r 7 ) 

kn 2 



(k,R) 



2mo 



75X1 E[M] (n-n)/f 15 + 



-^X2 E [kid] (2n? - - nf +2 )/f 15 + E {M} «i - n? +2 )/f aB + 



^Xs [(2 



£^5 



TV,) 7[ 12 + ([kj] N x - [k y 5] N y ) J 2 ri2 ] + 



a/6x6E{[W] ^ +2 }/f 15 + 75X7 E [{k<$} x N^/f 2 



APPENDIX B: MATRICES 

Matrices ff Q , used to obtain the H^^^ik, R) block of the r 8 ©r 7 two-band Hamil- 
tonian are: 



If 12 = 




-1 
1 




24 



Il 12 



I^ 15 = 



rri 5 _ 

y 



7^25 
V 






1 " 




n 

u 













J 


_ -1 


. 




—i 





n 












i 


1 





(J 


-V3 


-V3 








1 





-i2 




n 

U 











-il 





V3 








1 


-1 





- 


-V3. 










— i 




—i 



jr 2 5 





-2 
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